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Abstract: This paper reports the experimental and theoretical investigation 
of the Talbot effect beyond the paraxial limit at optical frequencies. Au hole 
array films with periodicity a 0 comparable to the wavelength of coherent 
illumination X were used to study the non-paraxial Talbot effect. Significant 
differences from the paraxial (classical) Talbot effect were observed. 
Depending on the ratio of a 0 1 A , the interference pattern in the direction 

perpendicular to the hole array was not necessarily periodic, and the self- 
image distances deviated from the paraxial Talbot distances. Defects within 
the hole array film or above the film were healed in the self-images as the 
light propagated from the surface. 
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The Talbot effect is a ubiquitous, lensless self-imaging characteristic of periodic arrays 
[l]encountered in a range of research areas from atom optics [2] plasmonics [3]. Under 
coherent illumination with wavelength X at normal incidence, an array with periodicity a 0 will 

form self-images at integer multiples of the Talbot distance z T = 2a„ IX [4]. In addition to 

self-images, patterns with periodicity of a 0 1 q can be observed at fractional Talbot distances 

pz T I 2q , where p and q are relative prime numbers [5]. Self-images and fractional images 

have been used in various applications such as lithography [6], interferometry [7], optical 
trapping [8], and array illumination [9]. At Talbot distances, self-images of the periodic array 
can be formed even in the presence of defects; this property has been used to generate defect- 
free self-images from a defective mask [10]. 

The self-imaging properties of the Talbot effect, however, are only valid for periodic 
objects under the paraxial approximation (a 0 IX »1) [11]. Previous work showed that lateral 
periodicity was not necessary for achieving longitudinal self-imaging, and certain non- 
periodic structures could generate self-images [12]. Beyond the paraxial limit, the periodicity 
in the x-y plane is not sufficient to ensure self-imaging in the longitudinal direction. Many 
applications based on the paraxial Talbot effect, therefore, cannot be extended to nanoscale 
structures with periodicities comparable to the wavelength of the incident light. Theoretical 
work has focused on the light patterns from ID gratings beyond the paraxial limit (the non- 
paraxial Talbot effect) [5, 13, 14], and some differences in the light patterns were predicted. 
For example, exact self-images only exist for certain ratios of a 0 1 X , and these self-images 
are expected to appear at positions smaller than Talbot distances. These predicted results, 
however, have not been confirmed experimentally because of challenges in fabricating large- 
scale periodic patterns that satisfy a 0 ~ X and difficulties in imaging the 3D light structure 

[15]. Recently, the non-paraxial Talbot effect was observed from surface plasmons launched 
by periodic groves on an Au film [16]; however, the fast damping of surface plasmons limited 
the study to only the first Talbot distance. 

Here we report an experimental and theoretical investigation of the Talbot effect beyond 
the paraxial limit. We measured and calculated the interference patterns formed by coherent 
light through hole arrays with a 0 on the same order as X. Different from the paraxial Talbot 
effect, the self-images were not periodic in the longitudinal direction. In addition, the 
measured self-image distance z R was different from the Talbot distance z T . Our theoretical 

analysis showed the deviation of z R from z T was related to the in-phase condition and strongly 
depended on the ratio a 0 1 X . The healing property of the Talbot effect was still observed in 
the non-paraxial regime. Because of constructive interference of plane waves transmitted 
through the periodic hole array, defect-free patterns could be observed at z R even in the 

presence of defects in the optical path. 

Scalar wave theory is traditionally used for the Talbot effect in the paraxial limit, where 
the light pattern through a periodic hole array at z = 0 (the initial structure) can be represented 
by a scalar wave if/ and Fourier expanded as a sum of plane waves in the direction of light 
propagation (z>0) [17] 

^.A#)=i;_i; MD ^(«)xexp[ic*+^*y+^d (i) 

r = (x, y, z) is the spatial position, X^- = XJ n is the effective wavelength in the medium with 
refractive index n, (a>) is the complex amplitude of the plane wave, and n\ and m 2 are 
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integers. The wave vector k = (fc™ 1 ,k^ ,&™"'" 2 ) = (2nn\ I a^^nn^ I a 0 ,k™ l,mi ) satisfies 

\k\ = 2x/A eff andk^" b >0. 

In the paraxial limit (a 0 »X e g), Taylor expansion of k z results in 

k""" 2 = (\k\ 2 - k? 2 - k^f 1 * 2n\\ - (m 2 + m 2 2 )A 2 ff / 2a 2 0 ] I A eff (2) 

If this approximation is applied to Eq. (1), the following relation can be derived 

y/(x, y,z + sz T , \ ff ) = exp[;2^iz r / X eff ]y/{x, y, z, \ ff ) (3) 

where z T = 2a Q I A eff is the Talbot distance and s is an integer. Thus, exact self-images will be 
observed at integer multiples of z T because the intensities satisfy 

I(x, y,z + sz T , l eg ) = I(x, y, z, l efs ) = |^(x, y, z, l 4f )| 2 
Beyond the paraxial limit, however, when a 0 1 X eff approaches unity, the Taylor expansion 

in Eq. (2) is not valid. The light patterns can still be calculated based on the plane wave 
summation, but the approximation about k z can no longer be used. To calculate the accurate 
phase change of the propagating plane waves (the exp[;7: z z] term in Eq. (1)), we consider k z 
independently for each (m^m^) plane wave beyond the paraxial limit. Since k z is imaginary 
for evanescent waves with k 2 x +k 2 y > (2k I A eff ) 2 , we therefore simplified the calculation of 

light patterns in the far field by neglecting the exponentially decaying evanescent waves. 
Thus, only the propagating plane waves, graphically depicted as points with coordinates 
(m 1 ,m 2 ) inside the circle with radius a 0 1 X eff in k space (Figs. 1(A)-1(B)), are considered. In 

addition, the condition for achieving exact self-images requires all the plane waves to be in 
phase at the same distance. If the center frequency (0,0,/c z °'°) is chosen as the reference, this 
condition can be written as: 

exp[/^zj = exp[/^°zj (4) 

for all propagating waves (m l ,m 2 ) , where z R is the self-image distance. We divided the non- 
paraxial Talbot effect into three different regimes based on the number of propagating waves: 
(i) a 0 1 X eff < 1 , (ii) 1 < a 0 1 \ ff < ^2 and (iii) a 0 1 \ ff > -J2 . 

In regime (i), only the center spatial frequency point in k space can propagate to the far 
field. Since there are no other plane waves inside the circle a 0 1 X eff to interfere with this center 

wave, no structured light pattern is observed. In regime (ii), four nearest neighbors are 
considered when calculating the light pattern in addition to the center spatial frequency 
(Fig. 1(A)). Because the four nearest neighbors have the same k z component, the self-image 
distance can be easily derived from Eq. (4): 
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Fig. 1. Self-image distances of the non-paraxial Talbot effect are different from the classical 
Talbot distances. (A-B)Fourier space spectrum of plane waves emitting from the hole array, 
where the circle with radius asJKg separates the propagating waves and evanescent waves. 
Light patterns from scalar wave calculations with (C) aJX e g = 1.230 (a 0 = 600 nm, Xgg = 488 
run), (D)ao/l e ff= 1.90 (a 0 =1.2 urn, l e ff = 633 nm). Experimental measurement with (E) atJKg = 
1.230, (F)ao/^g- = 1.90. Light patterns from FDTD simulations with (G) a<JX e g = 1.230, 
(H)aoM c# = 1.90. 

Equation (5) is similar to the self-image distances for ID gratings with l<a 0 1 X eff <2 
[13]. In regime (iii), additional plane waves with higher spatial frequencies contribute to the 
far-field light pattern (Fig. 1(B)). The increased ratio of a 0 1 X ef} results in an increase in the 

number of plane waves that need to be considered in Eq. (4). 

We simulated the light patterns at the different regimes above using the scalar wave theory 
(Figs. 1(C)-1(D)). Since transmission through the hole array is not considered in this model, 
we assumed that light was uniformly transmitted through the holes and not through the film. 
In experiment, however, the transmission is enhanced by surface plasmons and depends 
strongly on the material property and geometry of the structures (Figs. l(E)-l(F)).For a more 
accurate calculation and to compare with experiments, we carried out 3D finite difference 
time domain (FDTD) simulations (Lumerical® software package) .In the FDTD simulations 
(Figs. 1(G)-1(H)), we used plane wave illumination polarized along the x axis and a uniform 
mesh of 4 nm. Periodic boundaries were used in the x and v directions and perfectly matched 
layer (PML) was used in the z direction. The optical constants of Au were taken from Johnson 
and Christy [18]. Since we are only interested in the self-image distances, we only focused on 
the light patterns on a relative intensity scale; for all figures, yellow represents the strongest 
intensity. 

To compare the light patterns in the three regimes, we used Au square hole arrays with a Q 

= 600 nm and 1200 nm, and continuous wave laser illumination X e ff = 488 nm, 543 nm and 
633 nm. Large-area (>1 cm 2 ) Au hole array films with film thicknesses of 130 nm and hole 
diameters of 150 nm were fabricated using the PEEL technique [19]. Homogenous 
illumination was obtained by passing the laser light, delivered by an optical fiber, through a 
collimator. The 3D light patterns generated from the hole arrays were then imaged using a 
confocal microscope (Nikon® D-Eclipse CI) with z steps of 0.1 urn. The light was collected 
through a pinhole in the confocal microscope, which removed out of focus light and improved 
the resolution of the images. 
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In regime (i) ( a 0 1 A eff < 1 ), we tested the case a Q I X eff = 0.948 ( a 0 = 600 nm, X, 



633 



nm).As expected, no periodic structure was observed because there was no interference of 

plane waves (data not shown). In regime (ii) ( 1 < a 0 / <V2), we tested a 0 /A eff =1.230 

(Fig. 1(A)). The corresponding scalar wave simulation showed periodic self-images in the z 
direction (Fig. 1(C)) that agreed with experiment, where we measured the self-image 

distances to be 1 .2 urn, about 20% less than the Talbot distance z T = 2a^ I A eJf = 1 .48 urn 
(Fig. 1(E)). Based on the analytical expression given by Eq. (5), the self-image distance was 
calculated to be z R = 1 .2 um, which agreed with the measurement. The light pattern from the 
FDTD simulation also matched with our measurements (Fig. 1(G)). 

In regime (iii) ( a 0 1 A eff > -Jl ), self-images satisfying I(x, y, z + Zr) = I(x, y, z) can only be 

observed when the initial structure satisfies certain specific conditions [17]. Although exact 
self-images do not exist for the periodic hole arrays we used, we found that light patterns at 
certain z R do not deviate significantly from the initial structure. Since these differences are 
barely distinguishable from noise from experimental misalignment of the light source, we still 
designate these planes "self-image planes". For periodic hole arrays with a 0 1 X eff = 1 .9 

(Fig. 1(B)), we observed three self-images within 10 urn from the film surface in the 
measurement (Fig. 1(F)): the first self-image appeared at 4 urn, the second at 7.1 urn, and the 
third at 9.6 urn. Compared to the first and the second self-images, the third self-image was of 
weaker intensity, and deviations from the initial structure were also more pronounced. The 
measurement was very different from the paraxial Talbot effect, where the Talbot distance 

Z T = 2al I A eJf = 4.55 jum indicates two repeats within 10 urn from the surface of the film with 
the second repeat at z = 2z T = 9.\fim . We note that the measured self-images were not 

periodic in the longitudinal direction, and the fractional Talbot planes commonly observed in 
the paraxial Talbot effect were also missing. These differences can be explained by the 
breakdown of the paraxial approximation, which affects the interference patterns at both the 
Talbot planes and the fractional Talbot planes. Both scalar wave (Fig. 1(D)) and FDTD 
simulations (Fig. 1(H)) showed similar self-imaging behavior. 



1st nearest neighbor 
2nd nearest neighbor 
3rd nearest neigbor 
measurement 
Talbot distance 




"(Aeff 

Fig. 2. Difference between the measured first self-image distance zr and the Talbot distance Zt 
is related to ao/lgg-. The curves indicate the distances where the 1st, 2nd and 3rd nearest 
neighbors in k space are in phase with the center spatial frequency and approach the Talbot 
distance with increasing atJXeg. At the self-image distances (black squares), the deviations from 
the phase matching curves are relatively small. The error bars indicate the experimental errors 
in determining the self-image distances due to the noises. 
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To understand when the self-images distances z R approach the Talbot distances z T , we 
measured the distance of the first self-image for a range of different ratios a 0 1 X- and plotted 
its relative difference to z T as a function of a 0 1 /L. (Fig. 2). Our experiments showed that z R 
was always smaller than z T , which is consistent with the Taylor expansion theory that first 
order expansion of k z gives an upper limit estimation of the self-imaging distance. Since the 
self-images would form at distances where Eq. (4) was satisfied for all the plane waves 
simultaneously, we plotted out the solutions for each plane wave (m b m 2 ). As expected, exact 
self-images cannot be found because the distances where plane waves are in phase with the 
reference wave are different for different (mi, m 2 ). At the measured first self-image distances 
(black squares), deviations from the in-phase conditions caused the light patterns to be 
different from the initial pattern. These deviations, however, are small enough that the 
differences in the light pattern are barely distinguishable from the noises in the system. At the 
Talbot distance z T (green dotted line), the deviations are much more obvious; therefore, the 
light patterns at z T are significantly different from the initial structure. The same analysis can 
be applied to the second and third self-images. With increasing a 0 1 X eff , all the phase 
matching curves converged to z T ■ We estimate that when a Q I X eff >10, the non-paraxial 

Talbot effect will approach the paraxial Talbot effect. 

Besides resolving the 3D light pattern from a perfect initial structure, we investigated how 
the presence of defects would affect the transmitted light pattern. Understanding the impact of 
defects on the 3D light pattern can potentially benefit applications that rely on the Talbot 
effect, including lithography and imaging. First, we studied the influence of a 2D defect, a 
1.5-um diameter hole, within the a Q = 1.2um periodic Au hole array film in regime (iii) with 
a 0 1 X eff = 2.20. We oversaturated the intensity at the position of the defect so that the light 

pattern from the hole arrays and the defect could be measured at the same time. When focused 
to the hole array film (z = 0), the defect appeared as a bright spot approximately 1.5 um in 
diameter (Fig. 3(A)). The presence of the defect resulted in the disappearance of several 
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Fig. 3. 2D defects on the film were healed gradually in the self-image planes. Light patterns at 
self-image distances (A)z = 0 um, (B)z = 4.4 um, (C)z = 8.4 um and (D)z = 13.1 um. (E) the yz 
cross-section. As the light propagated in the z direction, the intensity of the defect decreased 
while the size of the defect increased. (a 0 =1.2 um, X e g = 543 nm). 

small bright spots corresponding to the holes in the hole array film. At the first self-image 
distance, some of the missing holes began to be restored (Fig. 3(B)). At the second (Fig. 3(C)) 
and third self-image planes (Fig. 3(D)), more missing holes were reconstructed, and the 
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intensity of the defect continuously diminished. This process where defects in the initial 
structure are healed in the self-images can be seen clearly in the y-z cross-section 
(Fig. 3(E)). In the presence of the 2D defect, the hole array still generated self-images, and the 
self-image distances were not affected. Effectively, the defect showed behavior similar to a 
point source emitting spherical waves, whereas the size of the defect became larger, the 
intensity decreased in the propagation direction. Light through the defect and the periodic 
hole array appeared to evolve independently from each other. At self-image planes farther 
than 10 um, the influence of the defect became very weak, and defect-free self-images were 
again observed. 

In addition to defects within the hole array film, we studied the influence of 3D defects 
away from the surface of the film. We used anisotropic Au pyramid particles [19] 
approximately 3 \xm in diameter and embedded them in a uniform poly(dimethylsiloxane) 
(PDMS) matrix as scattering defects. The pyramids were distributed randomly 1 1 urn away 
from the surface of the film. At the film surface (z = 0 um), we only observed the periodic 
hole array; the particle was not visible (Fig. 4(A)). When we focused to z = 11.2 um 
(Fig. 4(B)), two bright spots in the light pattern were missing because of scattering by the 
pyramidal particle. The missing bright spots were restored in the self-images as the light 
propagated in the z direction, and only one bright spot was missing at the self-image plane 2.3 
urn from the particle (Fig. 4(C)). As the propagation distance increased, the periodic pattern 
was completely restored as if there were no defect in the light path (Fig. 4(D)). In the y-z 
cross-section, the self-imaging property of the periodic hole array was unchanged, and the 
influence of the particle was limited to a very local region in the longitudinal direction 
(Fig. 4(E)). 




Fig. 4. Defect away from the surface of the film was healed in the self-images of the hole array 
pattern. The light patterns at (A) the film surface z = 0 um, (B) the position of the particles z = 
1 1.2 um(C) the self-image planes at z = 13.5 um and (D) the self-image planes at z = 29 um. 
(E) yz cross-section of the 3D light pattern showed the particle had a very local influence on 
the light pattern. (a a = 1.2 um, I = 543 nm, n = 1 .4) 



The healing effect observed for the two types of defects can be explained based on the 
lack of constructive interference of plane waves. Defects and the periodic hole array are 
composed of plane waves with different spatial frequencies, where the periodic pattern only 
has plane waves with spatial frequencies k = ( iTirn^ I a 0 ,2n:m 2 / a 0 ,k™"" h ).In contrast, finite- 
size defects can be decomposed into plane waves with many other different spatial 
frequencies. At the self-image distances, the hole array pattern reappears because plane waves 
with k = (2n:m 1 1 a 0 ,2n:m 2 1 a Q , k™ 1 ,mi ) are in phase and interference constructively. Most of the 
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plane waves composing the defects, however, will have large deviations from the in-phase 
condition and destructively interfere. With increased propagation distance, deviations from 
the in-phase conditions become larger, which lead to the gradual disappearance of the defects. 

In conclusion, we demonstrated the non-paraxial Talbot effect at optical frequencies. We 
found that self-images of the initial 2D periodic structure were not necessarily periodic in the 
longitudinal direction and that phase matching curves could be used to understand and predict 
changes in the non-paraxial self-image distances compared to the classical Talbot effect. 
Interestingly, the non-paraxial Talbot effect shared the same healing property as the paraxial 
Talbot effect. Defects with in the initial hole array structure as well as defects above the film 
surface were healed in the self-images of the hole array. Although this work was performed at 
optical frequencies, we anticipate our findings can be applied to other frequency ranges as 
well as other fields of physics. 
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